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Abstract 

The known solutions to the zero-sound dispersion equation are considered as placed on 
the physical and unphysical sheets in the complex plane of frequency. 

I. Introduction 

In this paper we present the known solutions to the zero-sound dispersion equation ob¬ 
tained in the Landau kinetic theory 0-0 and in the microscopic theory in RPA B-ffl 
in the complex plane of frequency. The solutions are calculated for the different values of 
coupling constant and considered on the physical and unphysical sheets. 

In recent years a large attention is paid to the problem of stability of the nuclear 
matter and the phase transitions in it 0. i. 0- To investigate the phase transitions, it 
is useful to understand what kind of excitations becomes amplified and leads to unstability 
of matter with the changing of the temperature or density or coupling constant. In the 
fermi-liquid the amplified solutions appear at F 0 < —1. When we take into account the 
meson exchange then the zero-sound modes undergo the drastic modifications and the 
damping solutions can turn into the amplified ones and lead to the appearance of the 
unstability of the description of nuclear matter. 

In this note the solutions to the zero-sound dispersion equation are presented. They 
are given in such a form that permits to trace the influence of the external field and meson 
exchange on these solutions later on. 

Following the papers 0.0 we write the dispersion equation in the form 

1 = FU 0R (uj, k). (1) 

In this equation T is the effective quasiparticle interactions, u and k are the frequency 
and the wave vector of the excitations. The operator II 0jR (u;, k) is the retarded zero order 
on the quasiparticle interaction self-energy part (polarization operator). There is a known 
analytical expression for this polarization operator 0-i 

It was shown in papers 0.0. that the energies of the zero-sound collective excitations 
are the poles of the two-quasiparticle Green function over the quasiparticle-quasihole 
channel. The two-quasiparticle Green function G ® is determined through the singlc- 
quasiparticle Green function G (1) and the vertex function T 
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G { 2) = G ( 1 ) G (1) + G ( 1 ) G ( 1 ) rG ( 1 ) G (1) . 


When we consider the interacting quasiparticles near the Fermi surface with the small 
momenta transmitted then the vertex function T can be given by 0 
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where a is a residue in the pole of CA 1 '. In this equation the constant interaction between 
the quasiparticles is supported. On the other hand, T is the two-particle scattering am¬ 
plitude that satisfy to the two-particle equation like Bethe-Salpeter one for the particle 
and hole [2]. 

The analytical singularities of the two-particle scattering amplitude in the complex 
plane of frequency are well known {3]. On the real axis of the physical sheet of u the 
stable excitations and cuts are placed . The cuts correspond to the excitation of the free 
particle-hole pairs. The damping excitations (resonances) are disposed on the unphysical 
sheets under the cuts. Besides, in the fermi-liquid at a strong attraction between quasi¬ 
particles F q < — 1 the imaginary solutions to the Eq. m appear. They are situated on 
the positive and negative imaginary axes of the physical sheet. These solutions break the 
causality conditions and point out the phase transition in nuclear matter. In the present 
paper we consider the stable, damping and increasing solutions of Eq.lfTj) for the different 
magnitudes of the quasiparticle interaction. 

Following to 0. 0 we determine the physical sheet of the complex plane of frequency 
on which the stable solutions of Eq. Q are situated. At some values of F the solutions 
are placed both on physical and unphysical sheets. These unphysical sheets belong to the 
frequency Riemann surfaces that are determined by the form of n 0R (cj, k). 

We take T as a constant scalar interaction F = F 0 Co. Here Co is a factor inversed to 
the level density on the Fermi surface for two kinds of nucleons, Co = tt 2 /(2mpp). Fo is 
a dimensionless coupling constant of the scalar quasiparticle interaction. 

We consider Eq. m in the kinetic theory of Landau 0-i and in the microscopic 
theory in the random phase approximation (RPA) j2]-jl]. The results of these theories 
coincide at the long wavelenths. For the fixed F 0 and Fermi momenta pp, the solutions 
to Eq. m are the linear functions of the frequency u on the wave vector k: uj — ks^ in 
the kinetic theory. In RPA there is a more composite dependence, and we have branches 
of solutions a ){k). 

In Sect.2 the analytical continuation of n oi? (u;, k ) on the unphysical sheets of the 
complex plane is realized. 

In Sect.3 the solutions of Eq. m for the attractive interaction Fo < 0 are presented. 
It is shown that the part of solutions to m is situated on the unphysical sheets (sheets / 
and I'). 

In Sect.4 the solutions for the repulsive interaction F 0 > 0 are presented. It is shown 
how after the overlapping of the frequency of the zero-sound excitations, oo s (k), and the 


2 



particle-hole continuum, the branch co s (k) goes to the unphysical sheet (sheets II and 

in 

In the paper we consider the symmetric nuclear matter at zero temperature and the 
equilibrium density p = p 0 , pp = 0.268 GeV. During the computations the mass of the 
quasiparticles is taken equal to m — 0.8 m 0 and m 0 = 0.94 GeV. 

II. The structure of polarization operator 

In this section the expression for a polarization operator is given. The formulae repeat 
the known expressions 1,0. But they are presented in a form that is more convenient 
for the analytical continuation of II 0R (u;, k ) in to from the physical on unphysical sheets. 
At the beginning we consider the expression for the polarization operator in RPA. 

Recall that there is a condition on the magnitude of the retarded operator Il n/? (u k ) 
in the complex plane of to j2j 


Il oiJ (—a;*) = (n 0H )*(a;). (3) 

The expression for 11° is taken from 0-i- Let us write it in the form that does not 
content the overlapping logarithmic cuts. The causal operator 11° can be presented as 
a sum of two terms. One term describes the excitation of the particle-hole pair in the 
medium and the second one corresponds to absorption of it 

n°(u;, k) — <j)(oo, k) + <f)(—co, k). (4) 

For the real to, the retarded polarization operator is determined by the following way 4] 

Il 0/? (u;, k) = (Re + i sign(uj)Im)U°(to , k). 


The expression for k) is written for the excitations in nuclear matter consisting of 
two sorts of nucleons with two spin projection. 

For the wave vectors 0 < k < 2 pp the term (j)(co, k) has a form 
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For k > 2 pp the term 0(u k) is a Migdal’s function Jjjj: 
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Consider the cuts of II 0 (a;, k ) in the complex plane of co. From © we see that there 
are two cuts for k < 2pp (denote them / and II). They are determined by the first and 
the second logarithms in © and are placed on the real axis at co equal to 


I : 0 < co 


kp F 
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It is easy to see that cp(co, k ) © is finite in the branch point com = kp F — \k 2 due to the 
cancelation of the infinite contributions of the first and the second logarithms. 


The cuts I and II can be considered as corresponding to the excitations of the different 
particle-hole pairs. The cut I describes the excitation of a hole with the energy and a 

particle on the Fermi surface. Then co = ~ and the points of the cut correspond 

to the change of the angle 6 between pp and k in the interval < cosO < 1. The cut 
II describes the excitation of a particle with the energy s^ F+ % and a hole on the Fermi 

surface. Then co = and the points of the cut correspond to 1 — ^ < cosd < 1. 

The cuts of the function 4>(—co, k) lie on the negative real axis symmetrically with 
respect to the cuts of (j)(co,k). Thus, II 0 (u;,fc) has four cuts in the complex plane of co 
which are shown in Fig. la. When u > 0 then the cuts stem from <p(u, k ) while for oj < 0 
they are caused by 0(— ui, k ). 

For the both of the logarithms in © the infinite-sheeted Riemann surface can be 
construct [Sj. The logarithms are the single-valued function in it. 


When k grows, approaching 2 p F , cut / becomes shorter and degenerates into a point at 
k = 2pp- For k > 2pp there are two cuts on the real axis in Eq. ©: III and symmetrical 
III'. 

///: TP£ + lt <a ,<tP£ + W (8 ) 
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Now we define the physical sheet according to the papers l> 0 In the long wavelcnths 
limit k —> 0, ^ —> const , n°(o;, k ) © has a form 
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We accept that on the physical sheet the logarithm in m has an imaginary part on the 
upper edge of the cut equal to — m JSj and, consequently, +iri on the lower edge. 

The corresponding retarded II 0 ^ is (we denote s = 


n °r? 


r jnpF 
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Substituting Eq. m in (HD we obtain the known form of the zero-sound dispersion 
equation in the kinetic theory 

1 + -W ^Ljn ( u, m + kpF \ = f ln l±I (12) 

F 0 2 kp F \<jjm — kpp J 2 s — 1 

In Eq. fjl 1 11 there is a cut on the real axis for (— ^ Looking at Eq.© we 

consider this cut as consisting of two cuts I and V (Fig.la). At the long wavelenths Eq.(|7jl 
comes to I : 0 < u < 

m 


III. Solutions at F 0 < 0 


In this section we consider solutions to Eqs. © and m at Fq < 0 in the complex 
plane of u. 


The family of solutions to Eqs. © and m at F 0 < 0 we denote by u s d■ The letters 
”sd” are connected with the words ’’spinodal decomposition”. It is known jj that the 
thermodynamic stability condition says that the matter is stable if the partial derivative 
of the pressure with respect to the volume is negative: v) T < 0- The P rocess which 

takes place in matter, when this condition is broken, is known as spinodal decomposition 
[7j. It is shown in [2, that there is a relation of the effective quasiparticle interaction to 
the partial derivative 


dP 

dV 


V 2 3m 


(1 + F 0 ). 


(N is the number of particles). The stability condition is broken at F 0 < —1. It looks 
reasonable to use the notation uj s d(k) for the solutions at F 0 < —1. It is shown in this 
section that all solutions at F () < 0 can be considered as belonging to the same family 
since they continuously turn one into another at the continuous change of F 0 . Therefore 
the notation u s d(k) is used for all solutions at F 0 < 0. 


III.A The definition of the unphysical sheets 

Let us turn to Eqs. © and © and define the unphysical sheets in the complex plane of 
u. Consider the cut I. We denote the first logarithm in © as ln(zi) = In (- ——-or ) • 

\Lom—kp F +J 

While the frequency u goes along the cut /©: u = (0, — — y~) then z 1 is changed in 
the interval Z\ = (0, — oo). The cut / in the complex plane of u corresponds to the cut 
along the negative real axis in the complex plane of z\ (Fig.lb). When we go to the cut 
/ from above (the arrow in Fig.la), this corresponds that we go to the cut from below in 
the complex plane of Z\ (the arrow in Fig.lb). Going over under the cut in Fig.lb we pass 
on to the unphysical sheet neighboring with the physical one on the Riemann surface of 
ln(zi). The magnitudes of ln(zi) on the neighboring sheet differ by a quantity of (— 2m) 
from the magnitudes on the physical sheet at the same Z\. 

Thus, going over under the cut / in Fig.la we pass on to the unphysical sheet of the 
complex plane of u that is placed under the physical sheet. Let us denote the unphysical 
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sheet with the same letter as a cut: /. The magnitude of In 


, P 

um—kpF~\- — 


on the sheet / 


differ from the magnitudes on the physical sheet by the quantity of (—27ri). As example, 


the value of In 


i p 

um-kpF + — 


on the lower edge of cut I is In 


value at the same point of u on the unphysical sheet is In 


i k 2 

um—kpF + — 


-;- w 

Lum-kpF-\- — 


T 7ii. But the 
- 7ri. Thus, we 


have the continuous changing of logarithm along the arrow in Fig.la. Then we conclude 
that H 0R in Eq.© changes continuously as well. The part of the unphysical sheet / is 
shown in Fig.lc by the shading with the right slope. 

An important remark is that moving on the frequency Riemann surface of one of 
logarithms we stand on the physical sheet for the other logarithms. In another words, the 
Riemann surfaces of logarithms in Eq.© are independent. 

Now we define an unphysical sheet related to the cut I' (denote the sheet by I'). The 
cut /' stems from the first logarithm in k ). We designate ln(z [) = —In 


= In ( “ m+k ^ 

l urn 


7 P 

-um—kpF-\- — 


|. Going over under the cut of ln(z[) we add the shift — 2m to the mag¬ 


nitude of ln(z[). Then, the logarithm on the sheet I' under the lower edge of the cut I' 


is equal to hi 


ujm+kpp — 4- 


+ 7 ri — 2ni. This expression is used for the analytical con¬ 


tinuation to the sheet I'. The part of the sheet I' is shown on the Fig.lc by the shading 
with the left slope. 

The values of II 0R (u, k ) on the unphysical sheets / and I' defined by such a way are 
connected by the relation 




(13) 


Thus, we made two unphysical sheets / and Let us show that the magnitudes 
of ir’^u;, k) Q on the negative imaginary axes of / and I' are real and are the same, 
according to m Then Eqs. (HD, © and their solutions coincide as well. 

On the imaginary axis we denote u = iuJi. On the sheet / the operator II 0R (a;, k ) 0J) 
has the following view on the negative imaginary axis 
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We leave only the Erst and the second terms in © because only these terms have the 
additional shifts on the sheets I and I'. This expression is the analytical continuation of 
ln(—Zi) — 7 ri to the imaginary axis. 
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On the negative imaginary axis of the sheet /' we have for n 0R (u;, k) 


n oiJ (iuy, k) = 


iujiin 

kp F 


mpF 


7T 


zujpn , 

1 - In 


nopn 


kpp icopn — kp F + i/c 2 


In 


iujjn + kp F — ifc 2 


- vrf + ... 


-tujpn 


and after the analogous transformation we obtain Eq. m as well. 


We can conclude that there are unphysical sheets / and I', that coincide on the 
imaginary axis. In the next section it is shown that the solutions of Eqs. ©, m at 
F 0 < 0 are disposed on imaginary axes both of the physical and unphysical sheets. 


III.B Solution in the kinetic theory of Landau 

At F 0 < — 1 the positive and negative symmetrical solutions to Eq. m are placed on 
the imaginary axis of the physical sheet (0. i). These solutions satisfy to the equation 
which can be deduced from Ea. (fT!Zl) (s = ij) 

1 + tt = 7 arctg- . (17) 

F 0 7 

The existence of the amplified solutions means that there is an unstability of the 
nuclear matter. 

Now we obtain the dispersion equation for the zero-sound excitations at — 1 < F 0 < 0. 
We consider Eq. m at long wavelenths and then substitute it into Eq. m- Then on the 
negative imaginary axis of the unphysical sheet we get an equation 

1 + 4- = ~1 ( ar ctg(l) ~ ^). (18) 

to z 

This equation was presented in pm where it was mentioned that there are no actual 
solutions of m- We really see that the solutions to Eq. m are placed on the unphysical 
sheet (/). 

In Fig.2 the solutions 7 of Eq. (ED at Fq < 0 are presented. The curves in this figure 
are the same as the curves in Fig.l in papers PI and [ 6 j (except the line ATI). The dashed 
domain means that the solutions are on the unphysical sheets. Looking at Fig.2 we see 
how the overdamping solutions on the unphysical sheet turn into the growing solutions 
and pass on to the physical sheet when the attraction between quasiparticles increases. 

Using the unphysical sheets one may obtain the symmetric solutions at — 1 < Fq < 0 
(the line AB). They are placed on the unphysical sheets (/ and F) that are above the 
physical sheet. In appendix A the construction of these sheets is presented. In Fig.2 the 
horizontal shading marks the sheet /. The dispersion equation ED on the unphysical 
positive imaginary axis has the same form on the sheets / and I'\ 

1 + 4- = —1 (arctg ( 7 ) + ^). (19) 

F 0 Z 
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This equation turns into Eq. m if to change the sign of 7 . Solutions of this equation give 
the curve AB in Fig.2. 


III.C Solutions in RPA 

In Fig.3 the solutions of Eq. (JTJ) obtained in the kinetic theory of Landau and in RPA are 
shown simultaneously at F 0 —~ 1.1, -1.2. The drawn solutions are on the upper semiplane 
of the physical sheet, i.e. they describe the unstable state of the nuclear matter. 

The kinetic theory gives the linear dependence u ~ k, that shows the unlimited 
undamping increasing of the frequency with k. In RPA the dependence is different: 
Imuj s d{k) reach the maximum and then decrease. At a certain value of k (which we 
denote by kf in ) the branch u> s d(k ) goes over under the cut to the unphysical sheet I. The 
values of the maximum (which is proportional to the growth rate) and kf in depend on 
parameters of matter F 0 , p F , m. 

In Fig.4 the branches u s d{k) obtained in RPA are presented. One sees that the 
branches change continuously with Fq passing from the physical to unphysical sheets 
I. At — 1 < F 0 < 0 c Vsd(k) are placed on the unphysical sheet completely. At F 0 < — 1 
the branches lie on the physical sheet at the momenta k within the interval k = (0, kfi n ). 
All solutions belong to the same family iv s d- 

Similarly to Fig.2 the symmetrical branch — u s d{k ) can be found. At k = kf in it goes 
over from the physical to unphysical sheets I. 

There are two questions related to Fig.4. It is shown how the branch go over the cut 
to unphysical sheet / at k — kfi n . But the cuts /, I', /, V exist when k < 2 p F . If a > s d(k) 
is on the unphysical sheet at k ~ 2 pp, where will it be at k > 2 pp when the cut is closed 
and teared off the physical sheet? One could show that uj s d(k) keep on to exist on the 
unphysical sheet. 

The second question is what happens with a J s d(k) when kf in > 2 p F I This takes place 
when Fq -C —1. Looking at Fig.4 we wait that at k — kfi n the imaginary branch u s d(k) 
go over to unphysical sheet but kji n > 2 pp and the cut / is closed. At kfi n > 2 p F there 
are two cuts III and III' (JHJ) on the real axis. They can be far from the point of origin 
if kf in is distinctly larger 2 p F . It can be shown that two imaginary branches u s d(k ) and 
—a } s d{k) turn to zero at k — kfi n and become real at k > kf in . Further, as k increase 
they move along the real axis in the different sides, reach the cuts III and III' and go 
over under the cuts. This behaviour of branches is a standard one and is describe, for 
example, in j<5j. 


IV. Solutions at Fq > 0 

The existence of zero-sound collective excitations at the repulsive quasiparticle inter¬ 
action Fq > 0 was predicted in the paper [I]. The solutions of Eq. 0 describing the 
zero-sound excitations in kinetic theory are presented in Fig.5. They are the same as in 
Fig.l in papers and ra. 


In RPA zero-sound excitations are described by the solutions of Eq. (HD• Usually they 
are denoted by u> s (k). The excitations propagate with a small damping till the overlapping 
of the frequency tv s (k) with the frequency of the free particle-hole pairs. The overlapping 
takes place at a certain wave vector k that we denote kd■ It is k — kd when the right 
point of the cut II meets a ) s (k d ) |3]. For k > k d the real solutions u s (k) that correspond 
to the stable excitations, turn into the complex ones. These complex solutions describe 
the damping excitations and are situated on the unphysical sheet II, Fig.6. 


The unphysical sheet II is made by the same way as in the previous section. The 
second logarithm in Eq.® has the value Ini—z 2 ) + m, where z 2 = k , PF , l ! 2 , on the 

OJTTL Kp p 2 K 

upper edge of the cut II and ln(—z 2 ) —ni on the lower edge. We define as the unphysical 
sheet II the sheet neighboring with the physical one where the magnitudes of ln(z 2 ) differ 
by the quantity of -\-2iri. As before going over from the physical sheet to the sheet II we 
have the continues changes of ln(z 2 ) and, consequently, of polarization operator. 

In Fig.6 the sheet II is marked by the checked shading. In this figure u s (k) is presented 
in the complex plane of uj. At the wave vectors k < kd, the real branch u s (k) is on the 
real axis. At k = kd, u J s (k) becomes damping, goes over the cut to the sheet II. 


In Fig.7 the different models of u s (k) are compared. The results are obtained at 
F 0 = 2, the corresponding wave vector is ^ = 0.51. The solid line stands for the kinetic 
theory. The dashed lines are for the real and imaginary parts of u s (k) in RPA. The 
imaginary part of excitations appears due to Landau damping :2]. 


The Landau damping can be calculated by the approximated equation [Ij. Let us 
denote a i r = Re(u s (k)) and uy = Im(u s (k)). Then at - C 1 the following approximate 
equation is valid 


u>i = sign(uj r )lmll 0 (uj r , k) 


(dRe{ n° 


l 


dcu 


( 20 ) 


The uji obtained with this equation is shown by the dotted line in Fig. 7. 


V. Discussion 


The effective quasiparticle interaction used in this paper is a simple constant interac¬ 
tion :5] 

T = C 0 (F + F'(tit 2 ) + G(aia 2 ) + G / (aia 2 )(rir 2 )), (21) 

where r, a are the isospin and spin Pauli matrices. The results presented above are the 
valid for all types of the interaction in Eq. (ED: scalar, isospin, spin and spin-isospin one. 

In the paper the unphysical sheets I and /' are defined. The family of the branches 
of solutions u> s d(k ) is placed (partly) on the imaginary unphysical axis where these sheets 
coincide. It seems that the remaining surface of sheets is not use. But we can turn to the 
consideration of the pion dispersion equation |5J in the complex plane of u. This equation 
is closely connected with the zero-sound dispersion equation in the spin-isospin channel. 
It is appeared that the solutions responsible for the pion condensation are situated on the 
sheets I and I'. When these solutions at a certain conditions go over to the physical sheet 
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(analogously a ) s d(k) in Fig.4) then the physical values become infinite. This is interpreted 
as a phase transition into another state containing the pion condensate in the ground 
state |5].[H]. 

The author thanks S.V. Tolokonnikov for the constructive critics and M.G.Ryskin for 
the fruitful discussions. The work is supported by the RFBR (a grant number 03-02- 
17724). 


Appendix A 


Here we construct the unphysical sheets / and /'. The solutions of Eq.(5) that are 
presented by the curve AB in Fig.2 are placed on these sheets. 

To pass on to the sheet I we do the analytical continuation in u> of ln(z) = In 


j z~2 

um—kpF + -^- 


from the lower edge of the cut / upwards, adding 2m to ln(z). The sheet / is situated 
above the physical sheet. The value of In 


———tt above the upper edge of the cut 

v ujm-kpF+-J 

I is ln(—z ) — 7 tz T 2m. This expression is continued on the sheet I. 

In the similar way we construct the sheet I' making the analytical continuation of 


In(z') = hi 


um+kpF —2“ 


^ from the lower edge of cut I' adding 2ni to logarithm. 
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1 Figure captions 


Fig.l. a)The complex plane of frequence. The cuts of fl 0 ^ Eq.Q are shown. Numbers 
stand for the cuts ©. b)The complex plane of z The cut corresponds to the cut / 
in Fig.a. c)The complex plane of frequence. The shading marks the unphysical sheets 
described in the text. 

Fig.2. Solutions u> s d(k ) to Ea. dT^l) at F 0 < 0. The variable 7 is 7 = ^/m uj sd . The 
shading with the right slope marks the unphysical sheet I. The horizontal shading marks 
the unphysical sheet I. 

Fig.3. The comparision of the solutions in the kinetic theory and in the RPA. The solid 
(dashed) lines correspond to solutions obtained at F 0 = —1.2 (F 0 = —1.1). 

Fig.4. The branches of solution uj s d{k ) obtained in RPA at different values of F 0 . The 
curve (1) corresponds to F 0 = —0.4; (2) F 0 = —0.9; (3) F 0 = —1.02; (4) F 0 = —1.1; (5) 
F 0 = —1.2. The shading marks the sheet /. 

Fig.5. The solutions to Eq. m u s at F 0 > 0. Two symmetric solutions are presented. 

Fig. 6 . The complex plane of frequence. The solutions oy(fc) to Eq.Q obtained in RPA 
are presented. The curve 1 is calculated at F 0 = 1 and the curve 2 at F 0 = 2 . The wave 
vector k d marks the point when the Landau damping starts at this F 0 : ^ = 0.13 when 
F 0 = 1, and ^ = 0.52 at F 0 = 2. 

Fig.7. The branch u) s {k ) at Fq = 2 is shown for the different models. The solid line is for 
the solutions of Eq. m- The real ay and the imaginary a* parts of the RPA solutions are 
presented by the dashed curves. The dotted line stands for the ay calculated by Eq. (HU). 
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